TOWARDS AUTOMATIC ANALYTIC EVALUATION OF 
DIAGRAMS WITH MASSES Q. 

Leo. V. Avdeev f\ J. Fleischer Q, M. Yu. Kalmykov ^ and M. N. Tentyukov Q 



Bogoliubov Laboratory of Theoretical Physics, Joint Institute for Nuclear 
Research, Dubna, Moscow Region 141980, Russian Federation 

Abstract 

A method to calculate two-loop self-energy diagrams of the Stan- 
dard Model is demonstrated. A direct physical application is the cal- 
culation of the two-loop electroweak contribution to the anomalous 
magnetic moment of the muon \{g — 2)^. Presently, we confine our- 
selves to a "toy" model with only u, 7 and a heavy neutral scalar parti- 
cle (Higgs) . The algorithm is implemented as a FORM-based program 
package. For generating and automatically evaluating any number of 
two-loop self-energy diagrams, a special C-program has been written. 
This program creates the initial FORM-expression for every diagram 
generated by QGRAF, executes the corresponding subroutines and 
sums up the final results. 
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1 Introduction 



Recent high precision experiments to verify the Standard Model of elec- 
troweak interactions require, on the side of the theory, high-precision cal- 
culations resulting in the evaluation of higher loop diagrams. For specific 
processes thousands of multiloop Feynman diagrams do contribute, and it 
turns out impossible to perform these calculations by hand. That makes the 
request for automatization a high-priority task. In this direction, several pro- 
gram packages are elaborated 0, || (|, || . It appears absolutely necessary that 
various groups produce their own solutions of handling this problem: various 
ways will be of different efficiency, have different domains of applicability, 
and last but not least, should eventually allow for completely independent 
checks of the final results. This point of view motivated us to seek our own 
way of automatic evaluation of Feynman diagrams. We have in mind only 
higher loop calculations (no multipoint functions). 

There exist several kinds of methods for evaluating multiloop Feynman 
diagrams with masses. The most direct and universal method is based on the 
Feynman parameter representation and subsequent numerical Monte-Carlo 
integration |6[]. The universality is its most essential advancement. Great 
progress has been achieved in implementing this method into combined au- 
tomatic systems like GRACE 0. But the convergence of the Monte-Carlo 
integration is rather slow, and there is no way to estimate the actual numeri- 
cal error. When the subintegral expression has kinematical singularities, the 
error may significantly exceed the estimate 0. 

In a semi-analytic method some integrations are done exactly in terms of 
special functions, and the dimension of numerical integration becomes lower 
[§], |J. The advantage of this method is the possibility to deal with tensor 
structures. However, an essential drawback is the difficulty of consistently 
performing renormalizations, because one has to stay in the integer dimen- 
sion. 



The method [10] of Taylor expansion of a diagram in external momenta, 
analytic continuation and Pade-like approximations allows one to recover 
the behaviour of the function in the whole complex plane of momentum 
variables. However, it would require the knowledge of rather many expansion 
coefficients in order to get a sufficiently precise estimate at the threshold. 

Therefore, we are going to use the asymptotic expansions of Feynman 
diagrams in small/large momenta and masses, where the orders of expansion 
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in a small ratio of parameters are completely collected even at the threshold. 
We are interested in the low-energy physics of the Standard Model. In this 
situation there exists a natural small parameter: the ratio of the character- 
istic scale of the process to the scale of the weak interaction, defined by M z . 
This provides a good convergence for asymptotic expansions of physical ob- 
servables and, in a number of cases, makes already the leading order sufficient 
for the existing precision of the experiments. On the other hand, one can use 
the dimensional renormalization, and then there are no problems with the R 
operation. 

We demonstrate here the functioning of a C program (TLAMM) for the 
evaluation of the two-loop anomalous magnetic moment (AMM) of the muon 
|(# — 2) . This piloting C program must read the diagrams generated by 
QGRAF |lTJ for a given physical process, generate the FORM [12| source 
code, start the FORM interpreter, read and sum up the results for the class 
of diagrams under consideration. Here, for the purpose of demonstration, 
we apply TLAMM to a closed subclass of diagrams of the Standard Model 
which we refer to as a "toy" model. 

Recent papers have reduced the theoretical uncertainty of the muon AMM 
by partially calculating the two-loop electroweak contributions Jl3|, |14|| . In 
some cases the following approximations were used: terms suppressed by 
(1 — 4sin 2 ©iy) were omitted; the fermion masses of the first two families 
were set to zero; diagrams with two or more scalar couplings to the muon, 
suppressed by the ratio were discarded; the Kobayashi-Maskawa matrix 
was assumed to be unity; the mass of the Higgs particle was assumed large 
as compared to Mz- 

All of these approximations, except possibly the last one, are well justi- 
fied and give rise to small corrections only. We consider it of great interest 
to study also the case ~ Mz- To perform this calculation is our main 
physical motivation. Apart from that, for technical reasons, it may be inter- 
esting to study the functioning of TLAMM by calculating all 1832 two-loop 
diagrams without any approximation. 

The calculation of the AMM of the muon reduces, after differentiation 
and contractions with projection operators, to diagrams of the propagator 
type with external momentum on the muon mass shell (for details see |T5|). 

The applicability of the asymptotic expansion |16|], |T?|] in the limit of 
large masses has to be investigated for diagrams evaluated on the muon mass 
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shell (i.e., p 2 = —rrvfy. Some diagrams already had a threshold at the muon 
mass shell before the expansion. In other diagrams, this threshold appears 
in some terms of the expansion. In dimensional regularization, threshold 
singularities (like any other infrared singularities if they are strong enough) 
manifest themselves as poles in e (in 4 — 2e dimensions). They ought to 
cancel for the total AMM. We check this for a closed subset of diagrams in 
the toy model. 

2 Large-mass expansion 

The asymptotic expansion in the limit of large masses is defined (|T7[ as 

F G (q, M, m, e) M ~°° £ F Gh (q, m, e) o T^F^, M, m\ e) (1) 

7 

where G is the original graph, 7's are subgraphs involved in the asymptotic 
expansion, G/j denotes shrinking 7 to a point; F 7 is the Feynman integral 
corresponding to 7; T Qi ^ is the Taylor operator expanding the integrand 
in small masses {m^} and small external momenta {g 7 } of the subgraph 7 
(before integration); "o" inserts the subgraph expansion in the numerator of 
the integrand Fq/ t The sum goes over all subgraphs 7 which (a) contain all 
lines with large masses, and (b) are one-particle irreducible relative to lines 
with light or zero masses. 

The following types of integrals occur in the asymptotic expansion of the 
muon AMM in the Standard Model: 

1. two- loop tadpole diagrams with various heavy masses on internal lines; 

2. two- loop self-energy diagrams, involving contributions from fermions 
lighter or heavier than the muon, with the external momentum on the 
muon mass shell; 

3. two- loop self-energy diagrams with two or three muon lines and the 
external momentum on the muon mass shell; 

4. various products of a one-loop self-energy diagram on shell and a one- 
loop tadpole with a heavy mass. 



4 



Almost all of these diagrams can be evaluated analytically using the pack- 
age SHELL2 For our calculation we have modified this package in the 
following way: 



1. There are no restrictions on the indices of the lines (powers of scalar 
denominators) . 

2. More recurrence relations are used, and the dependence on the space- 
time dimension is always explicitly reducible to powers of linear factors. 

3. A new algorithm for simplification of this rational fractions is imple- 
mented. These modifications essentially reduce the execution time (in 
some cases, down to the order of a hundredth). 

4. New programs for evaluating two-loop tadpole integrals with different 
masses are added. 

5. New programs were written for the asymptotic expansion of one- loop 
self-energy diagrams (relevant for renormalization) in the large-mass 
limit. 



3 The toy model 

As the first step, we concentrate on a "toy" model, a "slice" of the Standard 
Model, involving a light charged spinor ^, the photon A^, and a heavy neutral 
scalar field $. The scalar has triple (g) and quartic (A) self-interactions, and 
the Yukawa coupling (y) to the spinor. The Lagrangian of the toy model 
reads (in the Euclidean space-time) 

L = l d ^d^ + \M 2 :e:-^:e:-^:^:+\{d,A u -d u A,f 

+ ^-{d^f + (d + m)^ + zefi* - y$M (2) 

where e is the electric charge and a is a gauge fixing parameter. 
The main aims of the present investigation are the following: 
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1. Verification of the consistency of the large- mass asymptotic expansion 
with the external momentum on the mass shell of a small mass. In 
particular, we check the cancelation of all threshold singularities that 
appear in individual diagrams and manifest themselves as infrared poles 
in e. 

2. Estimation of the influence of a heavy neutral scalar particle on the 
AMM of the muon in the framework of the Standard Model. 

3. Verification of gauge independence (we use the covariant gauge with an 
arbitrary parameter a). 

In what follows we analyze in some detail the diagrams contributing to 
the AMM of the fermion in our toy model and specify the renormalization 
procedure. Apart from counterterms 40 diagrams contribute to the two-loop 
AMM of the fermion. After performing the Dirac and Lorentz algebra, all 
diagrams can be reduced to some set of scalar prototypes. A prototype defines 
the arrangement of massive and massless lines in a diagram. Individual 
integrals are specified by the powers of the scalar denominators, called indices 
of the lines. From the point of view of the asymptotic expansion method the 
topology of the diagram is essential. All diagrams of the toy model that 
contribute to the two-loop AMM can be classified in terms of 9 prototypes 
(we omit the pure QED diagrams). These prototypes and their corresponding 
subgraphs 7 involved in the asymptotic expansion, are given in Fig.|I|. In 
dimensional regularization, the last subgraphs vanish in cases 1,4,7, and 8, 
owing to massless tadpoles. 



4 Program TLAMM 

All diagrams are generated in symbolic form by means of QGRAF [II]. For 
automatic calculation we have created a special piloting program written in 
C. This program, called TLAMM, 

1. reads QGRAF output; 

2. creates a file containing the complete FORM program for calculating 
each diagram; 



6 



3. executes FORM; 



4. reads FORM output, picks out the result of the calculation, and builds 
the total sum of all diagrams in a single file which can be processed by 
FORM. 

The program has its own internal notation for topologies. In the problem 
g — 2 of a lepton in the Standard Model there are four different topologies 
(see Fig.|2|). Line number 1 is always assumed to correspond to a fermion line 
(a lepton or a neutrino); therefore, topologies b and c are distinguished. 
All diagrams are classified according to their local prototypes. The no- 
tation consists of a letter (the topology) and five (or four in B case) integer 
numbers specifying the masses on the lines 

for the zero mass: 7, v\ 

1 for a mass less than m M : e, u, d; 

2 for m M ; 

3 for an intermediate mass between m M and Mw'- s, r, c, b; 

4 for a heavy mass: W, Z,t, H . 

In the Feynman gauge a = 1 all pseudo-Higgs particles and the Faddeev- 
Popov ghosts are heavy except for one massless ghost. 

Each local prototype is calculated by means of the asymptotic expansion 
in heavy masses. For each diagram the corresponding FORM subroutine is 
called according to the local prototype. 

Identifiers for vertices and propagators and the explicit Feynman rules are 
red from separate files and then inserted into the FORM program. Because 
the number of identifiers needed for the calculation of all diagrams together 
may exceed the FORM capacity, the piloting program TLAMM retains for 
each diagram only those involved in its calculation. 

All initial settings are defined in a configuration file. The latter contains 
information on the file names, identifiers of topologies, the distribution of mo- 
menta, and the description of the model in terms of the notation that is some 
extension of QGRAF's. The program carries out the complete verification of 
all input files except the QGRAF output. 

There exist several options which allow one to process only the diagrams 

• explicitly listed by number; 

• of a given prototype; 
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• of a specified topology. 

There are also some debugging options. 

The asymptotic expansion of each prototype is performed by a separate 
FORM program. For efficiency of the algorithm the following points are 
essential: 

1. The result of the calculation is presented as a series in small parameters. 
Care is taken to avoid the production of unnecessarily high powers in 
intermediate results. 

2. For the evaluation of the Feynman integrals, it is necessary to reduce 
scalar products of momenta in the numerator to the square combina- 
tions which are present in the denominator. Most efficiently this is 
done by means of recurrence relations proposed by Tarasov [TjJ . 

As a demonstration of the functioning of TLAMM, the unrenormalized 
("bare") contributions of all the two-loop diagrams to the anomalous mag- 
netic moment of the fermion in the toy model are presented in Figs.|]-[7] to the 
leading order in m 2 /M 2 . In the presented expressions the factor of (27r) -4 is 
implied. During the calculations, each loop was divided by T(l + e) rather 
than multiplied by T(l — e) as is done in the MS definition. 

5 The results of the calculations 

Initially, we use the minimal subtraction scheme for the parameters of the 
Lagrangian (0). Afterwards, it is more convenient to re-express the running 
masses and M\ in terms of the physical pole masses m, M 2 of the particles 
(at the one- loop level, relevant here), 



m R {n 2 ) 
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(m 2 \ 2 f , M% /m 2 \ 3 /4 , M% 

+ 91 (l-^-iln^)l (4) 

thus terming to the on-shell renormalization of the masses. The scalar-field 
tadpoles do never contribute owing to the normal ordering of the Lagrangian. 
Anyway, their contributions would become of no consequence after everything 
is expressed in terms of the pole masses. As a result, the quartic interaction 
constant A falls out of the two-loop anomalous magnetic moment of the 
fermion. 

We have calculated the asymptotic expansion up to the 7th order in the 
ratio m 2 /M 2 and have convinced ourselves that all orders of the expansion 
are free from on-shell singularities and involve only the logarithms of the 
masses. For brevity, we present just two leading orders 
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and gn are renormalized (running) coupling constants. 

In quantum electrodynamics it is agreed to express the running charge 
of the electron efl(/z 2 ) in terms of the experimentally measurable physical 
charge e. The later is defined by the nonrelativistic Thompson limit of the 
Compton scattering, that is, as the product of the on-shell vertex function at 
the zero momentum of the photon [with the projection operator (—ip + m) 
In {~W + m ) ] by the on-shell wave-function renormalization constant (the 
residue of the fermion propagator at its pole). Both quantities are renormaliz- 
ation- and gauge-invariant but contain an infrared singularity which cancels 
in the product. Performing an analogous procedure for the charge of the 
muon in our toy model to the one-loop order we get 



The physical Yukawa charge could also be defined in terms of the on-shell 
Yukawa vertex. However, its evaluation is a difficult task in itself. On the 
other hand, in the Standard Model the Yukawa charge is usually related to 
the mass of the fermion rather than kept as an independent parameter. As 
a compromise, to define the physical Yukawa charge y, we use the following 
one-loop expression for the running charge: 



The calculation through the on-shell vertex would generally give a finite 
correction (independent of fi 2 ) to this formula. The running of the triple 
scalar interaction is inessential in the approximation that we consider. In 
terms of the physical charges, the anomalous magnetic moment (|5]) becomes 
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The Standard-model motivated values for the Yukawa and triple scalar 
interactions are 



1 
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Then, the estimate of the influence of a heavy neutral scalar particle on 
the anomalous magnetic moment of the muon is (besides the pure QED 
contribution) 
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This contribution is strongly suppressed and seems improbable to be exhib- 
ited in future experiments. 
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We conclude that 



1. the gauge independence of the two- loop contribution to the anomalous 
magnetic moment of the fermion has been verified; 

2. any threshold singularities in the on-shell asymptotic expansion of the 
diagrams contributing to the AMM do cancel; 

3. the correction due to the heavy neutral scalar is suppressed by the ratio 
of to heavy masses. 
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Figure 1: The prototypes and their subgraphs contributing to the large 
mass expansion. Thick, thin and dashed lines correspond to the heavy-mass, 
light-mass, and massless propagators, respectively. Dotted lines indicate the 
lines omitted in the subgraph 7. 
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Figure 3: The two-loop QED contributions to the AMM in the arbitrary 
gauge. 
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Figure 4: The two-loop AMM contributions proportional to y 4 . 
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Figure 5: The two-loop AMM contributions proportional to e 2 y 2 (to be 
continued). 
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Figure 6: The two-loop AMM contributions proportional to e 2 y 2 (contin- 
ued). 
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Figure 7: The two-loop AMM contributions involving the triple interaction 
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